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some given assigned arrangement. This result can be used for determining 
statistically whether stock judging is more than a matter of chance. 

O. S. Adams, Secretary-Treasurer. 



ON CERTAIN PROPERTIES OF MAKEHAM'S LAWS OF MORTALITY 

WITH APPLICATIONS. 1 

By H. L. RIETZ, University of Iowa. 

1. Introduction. Early in the development of the mathematics of life insur- 
ance, the calculations became laborious. This situation led to efforts to discover 
laws of mortality from which monetary values for life insurance could be calcu- 
lated. As early as 1725, De Moivre in his Annuities upon Lives, put forth his 
famous hypothesis that the decrements d x of the number living l x at age a; in a 
mortality table is constant for all values of x. Very little data were available to 
verify this hypothesis, but even the data of that time should have made it 
appear that this assumption could at most be a rough approximation. The 
De Moivre hypothesis is, of course, equivalent to a statement that the number 
living l x at age x can be given by a linear function of x. 

Just about one hundred years later in the Philosophical Transactions of 1825, 2 
Gompertz published the so called law of mortality which bears his name. Before 
discussing the hypothesis of Gompertz, it is desirable to explain the meaning of 
the expression "force of mortality." This expression is in common use by 
actuaries at the present time in connection with the doctrines of Gompertz and 
in other connections. It should perhaps be stated that Gompertz did not use 
the expression " force of mortality " in his papers, but used instead the expression 
" intensity of mortality " which has no technical meaning at present. 

The "force of mortality" ix x at any age x may be defined as the rate of change 
of the average number living l x of an indefinitely large class of persons of age in 
the neighborhood of x, per individual of the group l x . That is to say, if the 
function l x has a derivative with regard to x, 

1 dl x 

the negative sign being used so that ju x is positive. 

In other words, we conceive a continuous function l x which is proportional 
to the average number living at age x, out of an indefinitely large group of lives 
of age in the neighborhood of x. Then the force of mortality is defined as the 
derivative of this function with regard to age or time divided by l x . Obviously 

1 Read under slightly different title before the Mathematical Association of America, Sep- 
tember 6, 1920. 

2 "On the nature of the function expressive of the law of human mortality and a new mode 
of determining the value of life contingencies," pp. 513-583. 



1921.] makeham's laws of mortality. 159 

l x is a decreasing function. The negative sign is placed before the derivative 
dlx/dx so that n x is positive. Gompertz gave a rather wordy discussion of his 
hypothesis showing its reasonableness. 

Briefly, his hypothesis may be stated by saying that the rate of change of the 
force of mortality at any time is proportional to the force of mortality. That is, 

dftx 



or 



dx- Cfi * 



V-x = B(f, (1) 



where B and c are parameters. Then from (1) and the definition of force of 
mortality, we have 

-rr-Bf. (2) 

l x dx 
Integrating, we have in the usual form of Gompertz's law, 

l x = kg*, (3) 

where k, g, and c are the parameters. 

It is of some interest to know the view of Gompertz in regard to the real 
character of his derivation after he gave such a lengthy discussion of the reason- 
ableness of his hypotheses. He remarked that " this equation between the num- 
ber living and the age becomes deserving of attention not in consequence of its 
hypothetical deduction, which in fact is congruous with many natural effects, as 
for instance the exhaustion of the receiver of an air pump by strokes repeated at 
equal intervals of time, but it is deserving of attention because it appears cor- 
roborated during a long portion of life by experience." In the paper of 1825, 
Gompertz applied his function to portions of the Northampton Tables, Carlisle 
Tables, Desparcieux Tables, and certain Swedish Tables. 

In the Philosophical Magazine for 1839, 1 De Morgan showed that if the law 
of Gompertz applies to the whole of life, then Simpson's Rule, which has often 
been used as an approximation, to simplify joint life annuity calculations, would 
be rigidly correct. 

In 1860, Makeham suggested 2 his first modification of the law of Gompertz 
and in 1865 published a valuable paper setting forth the development of his 
suggestion. The hypothesis back of Makeham's first modification is that the 
force of mortality is given by 

iXx = A + B<f. (4) 

That is, the force of mortality consists of two parts — the one part increases in 
geometrical progression with age and the other is a constant. The addition of 
the constant seems reasonable when we recall that for certain diseases and 
accidents, the tendency to death seems practically independent of age. 

1 "On the rule for finding the value of an annuity on three lives," vol. 15, p. 337. 

2 Journal of the Institute of Actuaries, vol. 8, 1860, pp. 301-310; vol. 13, 1867, pp. 325-358. 
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Then from (4) and the definition of force of mortality, we have 

and by integration we have for Makeham's first modification 

l x = ks*g°*, (6) 

where k, s, g, and c are the parameters. 

The determination of these parameters to fit a given table of mortality under 
some criterion such as that of least squares is a problem of interest in attempts 
at accuracy and economy of time. This problem has received considerable 
attention in recent literature. It is perhaps sufficient for our purpose to state 
that the main methods used are : 

(a) Method of averaging results obtained from using four values of l x to corre- 

spond to four assigned values of x. 

(b) Method of moments. 

(c) Method of least squares. 

It may be of interest to write down the values of the parameters for some well- 
known table. For the American Experience Table, Arthur Hunter found the 
following: 

logio g = - .00013205 

logio c = .04579609 

logio s = - .003296862 

log M * = 5.03370116 

With respect to fitting tables well with the Makeham function l x = ks*g c ", 
there are many examples of tables that have been fitted well enough for practical 
purposes. As a striking example of closeness of fit, I wish to refer to a graduation 
of a rural life population table 1 for the registration area of the United States 
by C. H. Forsyth and the writer. The resulting curve crossed the curve of the 
given data eleven times in the range from age 19 to 85 which we treated. 

2. The advantages of a Makehamized table for joint life and survivorship in- 
surances. It is the main purpose of this paper to exhibit certain features of the 
application of Makeham's function in the treatment of joint life and survivorship 
insurances. It is fundamental in joint life and in survivorship insurances and 
annuities that we be able to find in the simplest way possible the value of joint 
life annuities to two or more persons of any assigned ages that may arise in the 
applications. We wish to show how it becomes possible with a Makehamized 
table to replace a cumbersome and laborious method of finding joint life annuities 
by a simple and elegant plan. 

1 Record of the American Institute of Actuaries, June, 1911, p. 19. 
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The laborious character of joint life calculations, when a table is not Make- 
hamized, should perhaps be indicated first. This can be shown with sufficient 
clearness by an illustration with only two lives of ages x and y. The joint life 
annuity for (a;) 1 and (y) is given by 

_ Vlx+lly+l + tflx+zly+i + l^lx+sh+S + ' ' ' ,-v 

Q'xy 7 7 , \* J> 

V X ly 

where a symbol I with a subscript means the number living in the mortality table 
at the age indicated by the subscript, and v is the discount factor. That is, 



v = 



1 + t' 
where i is the rate of interest. 

The formula (7) can be extended in an obvious manner to include three or- 
more lives. For only two lives, in De Morgan's forms, we introduce commutation! 
symbols as follows : 

Uxy == vxvy® , 

C X y = Qxly ~ lx+lly+l)v^^l\ 
— Vl)xy ~ Dx+i:y+i. 
Nxy ~ Dxy + Dx+1 : y+1 + D x +2 ; y+2 + * - " , 
M X y = Cxy + Cj;-)-! : y-|-i + C x +2 : y+2 + * • * , 
Sxy = Nxy + Nx+1 : y+1 + Nx+2 : y+2 + " * • , 
Rxy = M X y + Mx+i : y+i + M x +2 : y+2 + 

To exhibit the tabular values for any one of this set of symbols, say for D xy , 
for all combinations of ages from 10 to 95 inclusive would require a double entry 
table with (86) 2 tabular entries if the rates of mortality of the lives of (x) and {y) 
are given by different mortality tables, as in the case of a husband and wife 
when we use McClintock's Annuitants' Tables 2 (Male and female). If the rates 
of mortality of x and y are given by the same mortality table, the (86) 2 values 
would involve duplicates since D xy = D yx , and the number of necessary tabular 
entries would be 87^2. 

For three lives of (x), (y), and (z) with rates of mortality from different mor- 
tality tables, there would be required (86) 3 tabular entries to give D XV z- If the 
mortality rates for (x), (y), and (z) were given by a single mortality table, the 
number of necessary tabular entries would be 8 sC 3 . 

It becomes clear, without carrying the illustration to a larger number of lives, 
that, unless a very extensive joint life business were to be done, it would be 
impractical to prepare such tables. 

1 A symbol (x) is an abbreviation for "a person aged x." 

2 M. M. Dawson, Practical Lessons in Actuarial Science, New York, vol. 2, 1905, pp. 315-325. 



162 makeham's laws of mortality. [Apr., 

We should ordinarily find it more practical to treat each separate case by a 
method of approximation which consists in finding the approximate sum of the 
series of terms in (7). The method of approximation very often adopted con- 
sists in applying formulas derived from the Euler-Maclaurin formula of finite 
differences. Even with these methods of approximation, this plan is generally 
somewhat laborious. 

3. Lives of unequal ages can be replaced by lives of equal ages. We shall 
now attempt to give a brief exposition of the method that can be substituted for 
the rather laborious method of finding joint life annuities, when the table of 
mortality is a Makehamized table. 

In the usual notation, let 

be the probability that each of m persons (x{), (a; 2 ), (x m ) will live t years. Under 
Makeham's law, l x = hs x g c " for all values of a; in a certain interval. Then we 
have 

tP*. = 's l f m( - ct - l) . 
Then the value of the joint life annuity to (x{), (a; 2 ), • • •, (x m ) is 

Q>XiXf-X m = / V'P%iX 1 ,-»X m ~T~ V iPxiXf ■ ■%„ ~T~ " 3PxiX 2 ..-X m ~T~ ' '''T* tPx\Xz—X^T '} K°) 

= 2-i v tPxix s ...x m • • •> (") 

«=1 



z 



— y i S Bi,U(*i+t , *+'"+<*»)(t , -i) QO) 



where oi is the highest age in the mortality table. 

Let £ be the equal age of m persons that may,. for the purpose of finding the 
value of the annuity replace ages x\, x 2 , • • • , x m . 

Then t p% = s ( <7 c (c<_1) and the probability that m persons each of age £ sur- 
vive t years is 

(tP t ) m = s mt g mci< - ct - 1 ' 1 . (11) 

If we determine £ so that 

mc ( = d" 1 + c x ° + h <f, (12) 

we may clearly write in place of (10), 

<W-* m = £ vW" = s mt g mc (c ' _1) - (13) 



The solution of (12) for £ is generally obtained by the use of a table of forces of 
mortality. This plan is easily shown as follows : 
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From l x = k^g ", we have 

M * = ~~ T dx = ~~ log S ~ e log g log °' ^ 

From (12) and (14), we have that 

m^ = /**! + M* 2 + • • • + /**„• (15) 

Then the value of £ that satisfies (15) also satisfies (12). As tables of forces of 
mortality usually accompany a Makehamized mortality table, we solve for £ by 
interpolation from such tables. With m equal lives substituted for the m un- 
equal lives, there would be required for tabular entries of a commutation symbol, 
say for a D^ • • • , to cover all ages from some age such as 10 to an age such as 95 
inclusive, simply 86 values instead of the impractically large number of values 
mentioned above. 

4. Law of uniform seniority. For simplicity, let us consider first two lives 
(x{) and (x 2 ). The equivalent equal age £ is given by 

2c* = c* 1 + c x \ (16) 

If (x{) is the younger of the lives, then let x% = x\ + a, and we have 

2c* = e*(l + c a ). (17) 

From (17) we have 

£ = xi + k, 

where k depends only upon the difference of the ages a and not on the actual 
ages xi and x%. Hence, for a given difference between ages, we add the same 
number to one age to get the equivalent equal ages no matter what the ages xi 
and x-i are, 

To look at this in the general case of m lives we may note that if 

mc * = c* 1 + e* 2 + h c*™, (18) 

then 

mc l+t = c *i+* _|_ c *z+t _j_ 1_ c x„+t (19 ) 

for t years later, where t is any number. Thus, if we have obtained the equivalent 
equal age J at any date, then £ + t is the equivalent equal age for a date t years 
later. That is, the addition of a given number of years to each of the ages 
x\, x<i, • • • , x m will add the same number to £. This property of Gompertz's and 
Makeham's functions is known as uniform seniority. The applications of Make- 
ham's formula are much simplified because of this property. 

5. A single life may replace m lives by an appropriate change of the interest 
rate. For certain purposes, it is desirable to work with a single life $ for which the 
cost of an annuity would be the same as for joint lives (x\), (ar 2 ), • • • , (xm). When 
Gompertz's law holds, the replacement by a single life is possible without involving 
a change in the interest rate. But with Makeham's law, we proceed as follows. 



164 makeham's laws of mortality. CApr., 

We confine our attention to two lives of (x) and (y) since this appears simpler 
than m lives, although the plan can at once be extended to any number of lives. 
Recalling that 

tVx = 5 y*(«-i>, (20) 

and 

tVv = sY^- l) , (21) 

let £ be the single life to replace (x) and (y) in joint life probabilities. Then 

tp t = sy* (c<_1) - (22) 

From (20) and (21), 

tV*v = W+*x*-» = s l tPi, (23) 

where £ is given by the equation 

c £ = (f -f c v. (24) 

But the value of a life annuity is given by 

a X y = Vp xy + V 2 2 Px U + $ sPxy + •••+»' tPxy + ■ ' • (25) 

Let s» = v' be the modified discount factor. This indicates the change in 
interest rate to which we refer above. Then from (23), 

v'tPxv = rf'tPe (26) 

Hence 

a xy = ®'Pf + v' 2 2 p ( + v' 3 s p ( + •■•, (27) 

and we obtain the joint life annuities from the single age £. 

For special cases, and unusual rates of interest, the actuary is likely to find 
the substitution of a single life more convenient than the substitution of lives of 
equal age in number equal to the given number of lives. However, when com- 
plete tables are available, the plan that does not require a change of interest rate 
is likely to be preferred. 

6. Makeham's second modification of the law of Gompertz. A second modi- 
fication of the law of Gompertz was developed by Makeham in 1890. 1 

It may be recalled that the first modification assumed that the force of 
mortality 

Hx=A + B<?, 
where A is a constant. 

It is merely another step in the direction of generality to assume that fi x is 
the sum of the term Bc x used in the hypothesis of Gompertz, and a linear function 
A + Dx. We thus start from the assumption 

(x x = A + Dx + B<f, (28) 

or 

-jft-A + D'+B'. 

1 Journal of the Institute of Actuaries, vol. 28, pp. 152-159, pp. 185-192, pp. 316-331. 
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Then 

l x = fo-Js-wm-** = ks^Y. (29) 

l x+t = Jfes^'rW-oy". 

tVxy = ^'r 2 ^ 3 ^^ 2 '^^* 4 " ^^'-"- (30) 

Let £ be the age of equivalent equal lives, where £ is defined by 

2 C * = (f + c v , 

in case such a replacement by equal ages is possible. 
Then 

In the calculation of a joint life annuity, 

t — lM 

a xv = J^v'tPzv (32) 

«=i 

To determine a discount factor »' such that 

tftVxy = «'*#«, (33) 

we simply substitute in (33) from (30) and (31). This gives after some simplifica- 
tion 

or 

v' = OT 2 ^+")-«. (34) 

That is to say, using a discount factor v' given by (34) instead of v, enables 
us to replace lives of unequal ages by the same number of lives of equal age. 
Moreover, the law of uniform seniority holds for the same reasons given in other 
cases. 

This second modification has been used to graduate the new American- 
Canadian Mortality Table for Men, recently prepared at much labor and ex- 
pense, under the auspices of the Actuarial Society of America and the American 
Institute of Actuaries. 

It seems that further extensions of Makeham's functions expressive of human 
mortality could be made by using in place of A + Dx in (28) a polynomial of 
higher degree, but there would soon come a practical limit to the usefulness of 
such extensions. 

In conclusion, let me say that the main purpose in presenting this paper to 
the Mathematical Association is accomplished if I have succeeded in giving a 
general notion of those properties of Makeham's functions that make these 
functions of much practical value in the treatment of joint life and survivorship 
insurances and annuities. 



